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Solving Navier-Stokes equations coupled with a heat transfer equation 
using Bagarello’s approach and the Hankel transform 

Mahouton Norbert Hounkonnou, Villevo Adanhounme, and Jean Ghislain Compaore 


Abstract. In this paper, the dynamics of an incompressible fluid in a bounded connected 
domain, described by Navier-Stokes equations coupled with a heat transfer equation, is inves¬ 
tigated by a method inspired from the non-commutative strategy developed by Bagarello, (see 
Int. Jour, of Theoretical Physics^ 43, issue 12 (2004), p. 2371 - 2394). The solution of involved 
systems of partial differential equations is derived with the help of the unbounded self-adjoint 
densely defined Hamiltonian operator of the physical model and the Hankel transform. 
Keywords : Heat transfer equation, Navier-Stokes equations, Bagarello’s operator method, Han¬ 
kel transform. 
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1. Introduction 

Various models of motion are used to describe physical phenomena in fluid mechanics. See for 
instance the work based on Navier-Stokes equations reported in Refs. mu]. Unfortunately, due 
to the complexity of the involved systems of differential equations, most of these models are not 
exactly solvable. This explains the massive use of approximation and numerical schemes to obtain 
information on the flow phenomena such as the velocity distribution, flow pattern, pressure losses, 
and so on. To cite a few, see, e. g., m-M- 

In this paper, we show that the problem of the dynamics of an incompressible fluid in a 
bounded connected domain, described by Navier-Stokes equations coupled with a heat transfer 
equation, can be analytically solved combining both the operator method developed by Bagarello 
[1] and the Hankel transform. 

The paper is organized as follows. In section 2, the considered model is described. In section 
3, we obtain the exact solution to the heat transfer equation by the Bagarello’s non-commutative 
approach. The Navier-Stokes equations are solved by using the Hankel transform. Then follow 
the concluding remarks in section 4. 
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2. Mathematical formulation 


We consider cylindrical coordinates x = {r,(p,z) relative to the orthogonal basis (e^, 61 ^, 62 ) 
and incompressible viscous fluid flow with the 2 ;-axis symmetry. The fluid layer is confined between 
two parallel plates of horizontal lengths and Ly separated by a vertical {z) distance (d). The 
no-slip upper and lower plates are held at fixed temperatures Ti and T 2 respectively. A uniform 
volumetric heat flux H (with units power/volume ) is pumped into the layer. The governing 
equations of such a flow of an incompressible fluid in the standard Boussinesq approximation are 
Eglit[2] ,Lu[l] 


— -f u.Vu = —Vp + vV‘^\i-\-{gaT)\i in x (0, r), (1) 

dT 

— +U.VT = + j in (2) 

ot 

V.u = 0 in (3) 

with the initial and boundary conditions 

u|t=o = Uq, T\t-o = To in fl, (4) 

u = 0 in W = rx(0,T), (5) 

TU=o = Ti, T|2=d = r2, (6) 


where fl is a closed bounded set of with the Lipschitz continuous boundary F, u{t,x) = 
{u(t,r, z),0,v{t,r, z)) is the velocity of the fluid, p is the fluid pressure, T is the temperature, 
k = ( 0 , 0,1), ly is the viscosity, g is the acceleration of gravity along the z-axis (in the — k direction), 
a is the thermal expansion coefficient, k is the thermal diffusion coefficient. 7 = H/{pc)', p is the 
density and c is the specific heat capacity of the fluid. 


3. Analytical solution 

In this section, we investigate the analytical solutions to the above described model based on 
Navier-Stokes equations using an operator approach and the Hankel transform. 

As matter of commodity, we adopt cP /k as the unit of time, d as the unit of length and j k 
as the unit of temperature. The equations can then be put into the following dimensionless form: 

Pr-^(^+u.Vu) = -Vp-f V^u -I- RTk, (7) 

dT 

—+U.VT = V^T+l, (8) 

where Pr = vj k, is the Prandtl number and R = gaP^j(jpv) is the heat Rayleigh number. The 
boundary conditions now read: 

u = 0 in r =rx(0,r), T| 2 =o = T, T|2=i=0, (9) 

where F is the boundary obtained from F by replacing z = dhyz = l,T = {K/yd^){Ti — T2). 

We suppose that the solution of Cauchy problem exists and is unique. Then, we solve the 
problem by Bagarello’s approach [l], based on a suggestion drawn from quantum mechanics. It is 
based on the following: given a quantum mechanical system S and the related set of observables 
Os, i. e., the set of all the self-adjoint bounded (or more often unbounded) operators describing 
S, the evolution of any observable Y G Os satisfies the Heisenberg equation of motion (HOEM): 

^Y{t,x) =i[H,Y{t,x)]. (10) 

Here [A, H] = AB — BA is the commutator between A,Bg Os', H is assumed to be a densely 
defined self-adjoint Hamiltonian operator of the system acting on some Hilbert space T-L while 
the initial condition Y^ = h{x) is considered as an operator acting on the same Hilbert space TL. 
Then, the following statement holds: 
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Theorem 1. . A formal solution of the HOEM is Y{t,x) = where Y^ is the 

initial condition ofY(t,x) and H does not depend explicitly on time. Furthermore, if H is bounded, 

^ [it)^ 

we getY{t,x) = 

as : [A,B]o = B; [A,B]k = [A, [A,B]k-i], etc. 

d 


-[H,Y°]k, where [A,B]k is the multiple commutator defined recursively 


-Y{t,x) = 


dt 


gitHyO 


Proof. We have 
dt 

m 

Jt 

If H is bounded, Y{t,x) = can be expanded as 


d 




-Y{t,x) = 

-Y{t,x) = i[H,Y{t,x)\ 


dv 

e""“ Y^{—iH)e~ 




k=0 


□ 


□ 

Consider the finite Hankel transform of <i> defined by m 


h‘n')\t—0 



^Ht,0J3/2{l^n0df, 


Jo 


2 

a 


E 

n>l 


d3/2{t^nf) 




( 11 ) 

( 12 ) 

(13) 


where J 3/2 is the Bessel function of the first kind with the 3/2-order, pii, h-Si - ■ ■ t h-m ■ ■ ■ are the 
positive roots of the equation 

J 3 / 2 {tL£.o) = 0, fo= const. (14) 

Defining the stream function tf and the component of vorticity by 

\ dtp Idtp du dv 

^ r dz’ ^ r dr' ^ dz dr 
and taking the operator rot of the equation (O, we obtain [2] 


\dt dir, z) ) 

dT l ditP,T) 
dt r d{r, z) 

where 


5 f , diu;/r) \ 

(9r \ dr J 


3 a(a;/ 7 -) \l _ ^ 

dz\ dz J dr' 


1_9 / d^T 

r i9r \ dr ) ^ dz^ 


(15) 

(16) 


diA, B) _ dA dB dA dB 
d{r, z) ' dr dz dz dr ' 

Using the transformation f = ^/r'^ + z"^ and setting D(^,t) = uj{r,z,t)/r, T{^,t) 
the equations (ITSl) and (ITBl) become 


(17) 


Tir,z,t) - t, 


idn _ d^n Adn 

djff) ^ d^iff) 
dt df^ 


Rdf 


(18) 

( 19 ) 
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Using Bagarello’s approach, the solution of the equation m satisfying the initial condition can 
be expressed as: 




n>0 


n>0 


1 _ fn 

mt) = t+-Y,-A-[^To{r,z) 

^ n>0 


( 20 ) 

( 21 ) 


where A is the Laplacian operator. Using the change of variable fl = ^ 3/2$^ $ = in the 

equation (fTSl) we obtain 


Pr~^^ = 




52 $ 15$ (3/2)' 


dt oe e 

$it=o = 


dS. 


( 22 ) 

(23) 


where Uo(0 = (a;/r)|t=o- 

Applying the Hankel transformation to the equation (1221) and taking account of the boundary 
condition dO]), i. e. $(t,Co) = 0, we get 


Pr - Tr(t,M«), 

$|t=0 = ^(0,/in), 

_ rio 

Tit,fXn) = R 

Thus the solution of equation (l24l) satisfying the initial condition (l25l) can be expressed as 
¥(t, fin) = m - Pr f fin)dT 

Jo 

and 

[J'/,^(M!L^)]2'^3/2(/ine). 

Taking account of the expression of uj we obtain the equation with the unknown tp 
5V -2r%_3/2v- j 


(i$ 


and 


de 


(24) 

(25) 

(26) 

(27) 

(28) 

(29) 


and the derivative of stream function tp, components u and v can be defined by 

r.2 K 


dip -2r^ p - 3/2 ^ ^ 


. , - 2 r 2 /■« 

v[r,z,t) = -;r^ I 


n>l ^"^3/2^ 

-3/2 ^(tyd-n) j , 

d r7wrTAiI'^3/2(Mn??)d^. 


do [J2,/2^d-nio)Y 


Using Hankel transformation and Bagarello’approach, we provide analytical solutions to the non¬ 
linear Navier-Stokes equations and to heat transfer equation m and ([2]) satisfying the initial and 
boundary conditions ([3])-([5]). We can state the main result as: 
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Theorem 2. The velocity components u, v, the vorticity component to and the temperature 
T of the incompressible fluid flow described by the Navier-Stokes equations coupled with the heat 
transfer equation are defined by 


u{r,z,t) 


v{r,z,t) 


U}{r,z,t) 

T{r,z,t) 


-3/2 ^ 

1(11 ” 

2r 


^>1 ^'^ 3/21 


^'T-3/2 j , 

1 _ fn 

^ n>0 


(30) 

(31) 

(32) 

(33) 


4. Concluding remarks 

In this paper, we succeeded in solving the problem of fluid dynamics described by the Navier- 
Stokes equations coupled with the heat transfer equation. These equations were transformed in 
forms suitable to apply the Bagarello’s operator approach and the Hankel transform. Especially, 
we computed the velocity, the vorticity and the temperature of the fluid flow characterizing the 
fluid dynamics. 
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